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Abstract 

In this work we present non-linear models in two-dimensional space-time of 
two interacting scalar fields in the Lorentz and CPT violating scenarios. We 
discuss the soliton solutions for these models as well as the question of stability 
for them. This is done by generalizing a model recently published by Barreto 
and collaborators and also by getting new solutions for the model introduced by 
them. 
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Some years ago, Carrol et al [T] start to analyze the problem of Lorentz break- 
ing signature in field theoretical models. By now, there are a great number of works 
discussing a kind of symmetry breaking in many different physically interesting con- 
texts. For instance, in [2j it was discussed some impact over the standard model of 
this kind of symmetry breaking. Azatov ei a/ [3], in a recent work, have analyzed the 
spontaneous breaking of the four dimensional Lorentz invariance of the QED through 
a nonlinear vector potential constraint, Bezerra et al [1] have shown that a space-time 
with torsion interacting with a Maxwell field by means of a Chern-Simons like term can 
explain the optical activity in the sinchroton radiation emitted by cosmological distant 
radio sources, Lehnert et al [5j have verified the consequences over the Cerenkov effect 
of a Lorentz-violating vacuum and Bluhm j6] has made an estimative analysis of the 
Lorentz and CPT bounds attainable in Penning-trap experiments. In fact, along the 
last years a considerable effort has been drawn into this direction by many groups and 
in a variety of physical applications. On the other hand, the presence of topological 
solutions of nonlinear models is a matter of large interest and possible applications 
[3 [HI [9]. On the other hand, a natural place to apply these ideas is that of condensed 
matter non-relativistic ground, where the break of isotropy and homogeneity emerges 
quite naturally, due the material structure. 

As a consequence of the above arguments, it is natural to look for topological 
structures in CPT breaking scenarios. In fact, in a very recent work in this journal, 
Barreto et al [10] have introduced an approach capable of getting kinks in CPT violating 
scenarios. 

Here we are going to discuss a generalization of the work of reference [lOj, both 
by obtaining more general solutions for the models considered on that work and by 
generalizing Lorentz breaking Lagrangian densities. Particularly we obtain solutions 
which were absent in the reference [ID]. For this last accomplishment, we use a method 
recently introduced by one of us ^Ij. 

Models with Lorentz breaking terms usually leads to non-linear differential equa- 
tions, and one of the problems appearing as a consequence of this nonlinearity is that, 
in general, we loose the capability of getting the complete solutions. Here we extend 
an approach exposed in reference [11] which shows that for some two field systems 
in 1-1-1 dimensions, whose the second-order differential equations can be reduced to 
the solution of corresponding first-order equations (the so called Bolgomol'nyi-Prasad- 
Sommerfield (BPS) topological solitons [12j), one can obtain a differential equation 
relating the two coupled fields which, once solved, leads to the general orbit connect- 
ing the vacua of the model. In fact, the "trial and error" methods historically arose 
as a consequence of the intrinsic difficulty of getting general methods of solution for 
nonlinear differential equations. About two decades ago, Rajaraman [7] introduced an 
approach of this nature for the treatment of coupled relativistic scalar field theories 
in 1-1-1 dimensions. His procedure was model independent and could be used for the 
search of solutions in arbitrary coupled scalar models in 1-1-1 dimensions. However, 
the method is limited in terms of the generality of the solutions obtained and is con- 
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venient and profitable only for some particular, but important, cases [13]. Some years 
later, Bazeia and collaborators [14] applied the approach developed by Rajaraman to 
special cases where the solution of the nonlinear second-order differential equations 
are equivalent to the solution of corresponding first-order nonlinear coupled differen- 
tial equations. In this work we are going to present a procedure which is absolutely 
general when applied to Lorentz and CPT breaking systems, like those obtained from 
an extension of the ones described in [H] applied to nonbreaking versions appeared in 
^14j-pjj. Furthermore, we also show that many of these systems can be mapped into 
a first-order linear differential equation and, as a consequence, can be solved in order 
to get the general solution of the system. After that, we trace some comments about 
the consequences coming from these general solutions. 



1 BPS nonlinear Lorentz and CPT scenarios 

The two field model we shall study in 1 + 1 dimensions is described by the Lagrangian 
density 

^ = \{d,<Pf + \{d,xf - r (0, x)d,X - 9\<P. x)dA - V{<P, x), (1) 

where = 0, 1, x) and g^^cj), x) are vector functions with a prescribed functional 
dependence on the dynamical fields (f) and Xi and \^(0, x) is a potential term. 

Note that we can recover some usual Lorentz symmetry breaking models from ([1]) by 
choosing appropriately the vectors and g^. In particular, if /i = S2 and gi = Si x, 
one recovers the model introduced very recently by Barreto and collaborators (TU]. In 
fact, the first example we work out here is precisely this one, which we are going to 
show possesses an entire topological sector not considered in the work of reference |10j . 

If the potential V{(f), x) can be written in such a way that 

with W{(f), x) being any function of and Xy the energy density of the BPS states 
becomes 

^^^^=2i^ rf^ + ^^^'^'^V 2^ ^^^^^'^'^V ^ ^ 

with dW/dx = W^(j)' + W^x', where we have defined = = ^ and the 

prime stands for space derivative. 

From equation ([3]), we can see that the solutions of minimal energy are obtained 
from the following two coupled first order equations 

0' = W^{(f),x) - gi{<l),x) , 

(4) 

X = W^x(0,x)-/i(0,x) , 
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Finally the BPS energy is written, as usual, by 



EBPs = \W{(t),,Xj)-W{cP,,Xi) \ , (5) 

where 0j and Xi mean the i — th vacuum states of the model. Here, it is important 
to remark that the BPS solutions settle into vacuum states asymptotically. In other 
words, the vacuum states act as implicit boundary conditions of the BPS equations. 

It is interesting to notice that in the first order equations of motion (jlj) and in the 
energy density only the space components of the functional vectors and g^, fi 
and gi respectively, are present. 

From now on, in order to solve the equations (jlj), let us consider models for which 
we can write as a function of that is, In this situation, instead of applying 

the usual trial-orbit approach [14]-[19], we note that it is possible to write the following 
equation 

^'^ =dx= (6) 

where the differential element dx is a kind of invariant. In these cases one is lead to 

d4 ^ W^- gi 

dx W^-h' ^' 

Equation (jTl) is the generalization of the one studied in [11] to the case of nonlinear 
Lorentz and CPT breaking scenarios. It is, in general, a nonlinear differential equation 
relating the scalar fields of the model. If one is able to solve it completely for a given 
model, the function (p (x) can be used to eliminate one of the fields, so rendering the 
equations (jl]) uncoupled and equivalent to a single one. Finally, this uncoupled first- 
order nonlinear equation can be solved in general, even if numerically. 

We have found this method simpler than the method of the orbits broadly and 
successfully apphed to study the mapping of the soliton solutions and defect structures 
in problems involving the interaction two scalar fields. Despite of being simpler, the 
method applied here furnishes not only the same orbits than those obtained by using 
the method of the orbits appearing in the references [Il]-[in], but also some new ones 
as can be seen through a comparison with reference [11]. In the example worked out 
below one can verify that, this time, the mapping constructed here furnishes the very 
same orbits obtained in the reference [10]. Notwithstanding, we were able to find new 
solitonic configurations, not observed by Barreto and collaborators. 



2 The example of linear Lorentz and CPT breaking 

In this section we consider the particular model introduced in the work of Barreto et al 
[TU] in order to apply the method discussed in the previous section. In fact, we show in 
this example that the equation ([7]) can be mapped into a linear differential equation, 
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from which it is possible to obtain the general solutions for the soliton fields. In the 
case on the screen, the superpotential [TU] is written as 

W{<P,x)= 0-^0'-r0x', (8) 

and the Lorentz symmetry breaking terms in the lagrangian density ([1]) are chosen to 
be given by /i(0, = S2(p and gi{(p,C,) = Si x, such that equation ([7j) is rewritten as 



(9) 



and 



the equation Qj becomes 



# ^ (0^ -1) +r + SiX 
dx 2r(j)x + S2 (p 

where Si and S2 are constants. 

At this point one can verify that, performing the transformations 

x = C-|^, (10) 

b' = p+l + ^{2s^-S2), (11) 
dQ r Q r 

which is a typical inhomogeneous linear differential equation ^llj . The general solutions 
for the orbit equation are then easily obtained, giving 

02 - 1 = coC' + C - -^C + k for r ^ 1 and r 7^ ^, (13) 

2r — 1 r — 1 2 

0^ - 1 = -6 C In (C) + ci C + + k, for r = 1 (14) 

and 

02 - 1 = In (C) + K + C2 C' + fc, for r = ^, (15) 

where k = |^(2si — S2), b = S2 — Si and Cq, Ci and C2 are arbitrary integration 
constants. 

In general it is not possible to solve x in terms of from the above solutions, but 
the contrary is always granted. Here, with the aid of ([8]) and ffTOl) . we shall substitute 
the expressions of (x) obtained from f|T3l) . f|T^ and f|T5|) in the second equation (jl]), 
obtaining respectively: 



^ = ±2rcJl + CoC^ + -^C--^C + k, forr^l, r ^ 1, (16) 
dx V 2r — Ir — 1 2 
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dx 



±2r CVl-Kln(C) + CiC + e + A;, forr = 1, 

(17) 



^ = ±2 rCv^l + C' ln(C) + K + C2 C' + ^ , forr = ^. 



Barreto and collaborators [TO] have limited themselves to the orbits in which r 7^ 1 
and r 7^ 1/2 and the arbitrary constant Cq equals to zero or infinity. In the particular 
case with cq = they have found a lump-like profile for the field x(^) ^ kink-like 
profile for the field 0(x). By integrating the equation ( |T6l) and substituting its solutions 
into the equation (fTOl) we get the following solutions for the field x(x) 

y^^fx) = - fl8) 

^ (v^eT2v^^(^-^«) + C)2 - 2r' ^ ^ 

where xq is a constant of integration, A = 1 — b'^/Ar, B = r/(2r — 1), C = 6/(r — 1) 
and we have taken si = 0. On its turn the solutions for the field (pix) are obtained by 
substituting the classical solutions of the equation f|T6l) into the equation f|T3l) . namely 



The above solutions are valid if the parameters satisfy the conditions A > and 
7^ 4y4i?. The behavior of the above solutions are plotted in the figure 1 for the 
parameters r = 0.4 and b = 0.6. One can observe that in both pairs of solutions, 
(0+,X+) and ((/)_, X-), the field x{^) exhibits a lump-like profile and the field 0(x) 
a kink-like profile. This behavior is also found in many systems of two interacting 
solitons reported in the literature. 

More recently [TTj it has been shown that many models of two interacting solitons, 
very similar to this one with explicit Lorentz symmetry breaking that we are presenting 
here, can also exhibit kink-like behavior for both of the soliton fields, depending on the 
range of the parameters of the model. Inspired on this achievement, we have shown that 
it is also possible to have kink-like profiles for both of the fields, for particular values of 
the parameters r and b, in the model treated here. In fact if one takes 6 = 2(r — 1)/ y/r, 
which corresponds to one of the solutions with = 4AB, and r > 1/2 in the equations 
(fTSl) and ([T9l) we obtain the following forms for the fields 

R/ N 4(2r - 1) r - 1 
X± a; = \ ^ ^ ^, 20 



and 



^B^^. ^ ± 4(2r-l) .21) 
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In the figure 2 we present the behavior of the above kink solutions for r = 2. 

One could interpret these solutions as representing two kinds of torsion in a chain, 
represented through an orthogonal set of coordinates and x- So that, in the plane 
the type-y4 kink corresponds to a complete torsion and the type-i? corresponds 
to a half torsion, similarly to what has been done in [TT] . 

It is worth mentioning that the pairs of type-i? solutions have a BPS energy lower 
than that associated to the type-A soliton solutions. This can be shown by substituting 
the asymptotic values of the solutions in the equation ([5]), that is, for the type-A 
solutions we find E^pg = ^A\/~A, and E^pg = ^A\/A for the type-i? solutions. 



3 Generalized models 

In what follows, we will study a more general model contemplating a number of partic- 
ular cases which have been studied in the literature, including the previous and some 
other new ones. For this, we begin by defining the superpotential 



W^(0,x) = f + , (22) 



= ^ 

such that the equation (JTj) is given by 

(23) 



dx ixN ct>^x-fM^X) 

where = dF/d(j). The space-component of the functionals terms responsible for 
breaking the Lorentz symmetry explicitly, namely, fi{4>,x) ^ind gi{(f),x)y ^ire to be 
chosen more general than those of the model discussed previously and conveniently 
such that the integration of the equation fl23l) be possible. Based on the succesfuU 
generalization of models of interacting solitons also carried out in the reference [11] 
and in the development of the model of the previous section, a possible generalized 
model can be constructed by choosing 

and the following forms for the functionals fi{4>, x) and gi{(f), x), 



/i(0,x) = b^^x, 

^?i(0,x) = a0^-^x, (25) 

where is a positive integer number, A and 7 are constants and the parameters a and 
b can be thought as space-components of two-vectors pointing out in some preferred 
direction in space-time and the responsible for breaking the Lorentz symmetry. 
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The corresponding equation for the dependence of the field as a function of the 
field X is now given by 



dx 2 ■ ^ ^ 

Now, by performing the transformations 

^ = :^ + + 7 , (27) 

and 

^ = fiX-b, (28) 

we get 

dcr Act A^A A , , , , 

The above equation is very similar to the equation (fT2|) and can be easily integrated 
out. Its general solution in the case A 7^ /i and A 7^ 2/i is 

,, X{Nb-a) NX 2 A 

-(^^-7(;7rAy^+ 2,M2.-A) ^ 

where c is an arbitrary integration constant. The solutions for the equation fl29l) in the 
cases A = yU and A = 2/i can also be obtained, but we will not deal with them here. 

We substitute the equations fl27|) . fl28|) and fl30l) in one of the equations (jlj) to obtain 
the following first-order equation of motion for the field <j 



dx 



N 



2/i - A 



.2 I 2(iV& 
/i- 



A 



-1 N/2 



Nb'^ - 7/i + 2a6 



(31) 



This last equation can be solved analytically or numerically, depending on the values 
of the parameters. For the particular case with N = 2, 2b = a and c = we obtain 
very simple kink solutions for both of the fields (pix) and xi^)^ can be verified from 
the behavior of the solution for the field ^(x) 



^(x) 



Be 



Bx 



(7/i — 26^)/A > 0, and by substituting 



(32) 



in the 



where A = 2/{2fi - A) > and 5 
equations (|2H|). (I30l) and (1271). 

The construction of an even more general model which includes non-linear depen- 
dence on the field x{^) can be carried out by following the generalization proposed 



8 



in the reference [TT]. This can be accomphshed by choosing the following form of the 
superpotential 

W^^M(0,x) = ^0^X*' + i^(0), (33) 

where F{(j)) is given by the equation ([23]) and M is a positive integer. In order to include 
the terms responsible for breaking the Lorentz symmetry and to obtain a solution for 
the differential equation ([7]) it is reasonable to choose the functionals /i(0, x) 
gi{(f), x) the following forms 



/i(0,x) = &0V'-^ 

g,{<P,x) = a0^-V'- (34) 
With this generalization the equation ([7]) can be written in the form 

^ _ ^ ^-M ^ (/iiV - 2aM) 

dx {fi-b)X^ M{fi-b)X ^ ' ^ ' 

where ip = + 7. 

The equation fl5S]) is similar to the one which appears in reference [TT]. It admits 
the solution 



exp 



A(/i-6)M-2 



X 



(2-M) 



2M/(M~2) nN-2aM 2, 

5i + ..... ^^^.w — —X ' 



X 



(2-M) \2/(A^-2) 



M(M-2) 2\{^-h) \[M-2) 



X 



X 



(2-M) 



(M-2)' A(/i-6) (M-2) 



(36) 



where ci is an arbitrary integration constant and r(a, z) = t"' *dt is the incom- 
plete Gamma function. 



4 Conclusions 

We have been able to generalize a model presented recently in the reference ^Oj which 
incorporates the phenomena of solitons interactions and the Lorentz symmetry break- 
ing. The generalization has been carried out in two ways. We have found non-trivial 
classical solutions which exhibit kink-like behavior for both of the interacting fields 
and, consequently, with BPS energy lower than that associated with the usual so- 
lutions presented previously for the same model. Another interesting aspect of the 
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kink-like solutions rest on the study of the stability of the solutions against small time- 
dependent linear perturbation. At least for some models with only one scalar field, it 
has been shown in the reference ^Oj that models with kink-like solutions possess the 
stability of these solutions, on the other hand, models with lump-like classical solutions 
are unstable. For two interacting scalar fields the problem is cumbersome, even though 
the authors of the reference [10] have been able to show, based on very elegant and 
general arguments, that the solutions found there, even with lump-like configurations 
for one of the fields, are stable. We understand that the analysis of the stability car- 
ried out in reference [10] is valid for reference systems in which Bq = 0, where bo is the 
time-component of the two- vector responsible for the Lorentz symmetry breaking. For 
reference systems in which bo ^ the analysis has not been done. 

We have also proposed generalizations of the model of the reference [TO] by intro- 
ducing non-linear terms that break the Lorentz symmetry. This last generalization was 
possible thanks to the successful generalization carried out in the reference [H] which 
deals with a Lorentz symmetric two-dimensional model of interacting scalar fields. 
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Figure 1: Typical typc-A kink profile (for r = 0.6, 
to the field x+{^) the thick line to the field 0+ 



X) 



-- 0.4). The thin line corresponds 
Both were calculated for cq — 0. 




Figure 2: Typical type-S kink profile (for r — 2). The thin hue corresponds to the 
field x+(x) and the thick hue to the field 0+(x). Both were calculated for cq = 0. 
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